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ABSTRACT 

We report the discovery of a high C0( J=3-2)/CO( J=l-0) ratio gas with an arc-hke distribution 
( "high-ratio gas arc" ) surrounding the central star cluster of the supergiant Hii region NGC 604 in the 
nearby spiral galaxy M 33, based on multi- J CO observations of a 5' x 5' region of NGC 604 conducted 
using the ASTE 10-m and NRO 45-m telescopes. The discovered "high-ratio gas arc" extends to the 
south-east to north-west direction with a size of ~ 200 pc. The western part of the high-ratio gas 
arc closely coincides well with the shells of the Hii regions traced by Ha and radio continuum peaks. 
The C0( J=3-2)/CO(J=l-0) ratio, i?3_2/i_0j ranges between 0.3 and 1.2 in the observed region, and 
the i?3_2/i_o values of the high-ratio gas arc are around or higher than unity, indicating very warm 
(Ticin > 60 K) and dense (riHa ^ 10^~''cm~'^) conditions of the high-ratio gas arc. We suggest that the 
dense gas formation and second-generation star formation occur in the surrounding gas compressed 
by the stellar wind and/or supernova of the first-generation stars of NGC 604, i.e., the central star 
cluster of NGC 604. 

Subject headings: galaxies: individual (M 33) — ISM: Hii regions — stars: formation — ISM: indi- 
vidual (NGC 604) ~ ISM: molecules 



1. INTRODUCTION 

Giant or supergiant Hii regions (hereafter referred to 
as GHRs) are one of the most prominent objects in star- 
forming galaxies at the optical wavelength. Their sizes 
often reach a scale of a few 100 pc, and their Ha lumi- 
nosities are typically of the order of a few lO^^"**^ erg 
s~^, which corres ponds to a few 10 to a few 100 05 stars 
(jKennicuttI 1 1984f ) . The structure of GHRs is character- 
ized by (1) the presence of a central young star clus- 
ter and (2) extended shells and/or arc-like filaments sur- 
rounding the central star cluster. For example, 30 Dor in 
Large Magellanic Cloud (LMC), which is the most lumi- 
nous Hii region in the Local group, houses the compact 
cluster R136, which is cla ssified as a super star cluster 
(SSC) (|Hunter et al.l[l995l) . On the other hand, NGC 604 
resides in M 33, which is the second-most luminous su- 
pergiant Hl^£egwn_aftei_30 Dor, hosts a scaled OB asso - 
ciation (jHunter et al.lll996l : iMaiz-Apellaniz et al.ll2004D . 
These central star clusters appear to have formed during 
the initial stages of the formation of GHRs, and are ex- 
pected to have a strong impact on their natal molecular 
clouds due to their strong UV radiation, stellar wind, and 
supernova explosion. Therefore, GHRs provide us with 
an ideal environment to understand the clustered OB 
star formation process, and their impact on the ambient 
interstellar medium (ISM). These physical processes are 
also crucial in the evolution of starburst in galaxies. 

In this study, we present the ^^CO(J=3-2) and 
12CO(J=1-0) observations of NGC 604 in the near- 
est face-on spiral galaxy M 33 using the Atacama Sub- 
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millimeter Telescope Experiment (ASTE) 10-m and the 
Nobeyama Radio Observatory (NRO) 45-m telescopes. 
The structure of NGC 604 is complicated and contain- 
ing many shells, filaments, and arc-like structures; fur- 
thermore, there is a central star cluster surroun ded by 
arc-like Hii regions (iGomez de Castro et a l. 2000) and ^ 
200 massive OB stars (jConzalez Delgado &: Perez 2000 ) . 
Several radio components exist in the arc-like Hll regions, 
which are photoionized from the inside by the obscured 
massive stars, embedded in the further extended halo 
(jChurchwell fc Gosslll999f ). The atomic and continuum 
emission of NGC 604 in the far - infrared waveler igths is 
also bright (jHigdon et al.l 120031 : iHippelein et all 12003'). 
Its proximity (D = 0.84 Mpc; Freedman et al. 1991) 
and the favorable inclination angle of the host galaxy {i 
= 52°; Corbelh & Salucci 2000) allow us to clarify the 
detailed distributions of young stars and the associated 
molecular clouds within this extreme star forming region. 

Our objectives in conducting these molecular line ob- 
servations are (1) to investigate the spatial variation 
of the physical properties of molecular gas and (2) to 
understand their relationship with the GHR formation 
processes. Several researchers have previously reported 
on observations of M 33 in C0( J=l-0) emission, which 
can be coUisionally excited even in low-density molecu- 
lar gases such as n(H2) ^ 10^ cm~^, and a large molec- 
ular clou d complex associated with NG C 604 has been 
detected (IWilson fc S covillc 1992; Viallef ond et al.lll992t 
lEngargiola et al.l l2003i : iHever et al.i i200Ik Wilson et 
al. (1997) observed ^^C0{J=2-l), i3CO(J=2-l), and 
^^CO(J=3-2) emission lines, probing the denser com- 
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Minamimaki, Minamisaku, ponents of molecular gas (n(H2) > lO^"'^ cm'^), to 

ward two points in the NGC 604 region; they found that 
the kinetic temperatures of the molecular clouds asso- 
ciated with the Hii regions were systematically higher 
than those not associated with the Hii regions. Never- 
theless, no map has been created for these high-J CO 
lines thus far, and therefore the spatial variation of the 
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physical conditions of the molecular clouds and their re- 
lationship with the central star cluster and the surround- 
ing Ha shells/arc-like structures have not yet been well 
understood. 

Our new observations conducted using the ASTE 10- 
m and NRO 45-m telescopes provide high sensitivity 
i2cO(J=3-2) and ^'^CO{J=1-0) images of a 5' x 5' or 
1.3 X 1.3 kpc region located in the center of NGC 604 
with an effective spatial resolution of 25" or 100 pc. Us- 
ing these new images, we report the discovery of a high 
C0( J=3-2)/CO( J=1~0) ratio gas with an arc-like distri- 
bution (high-ratio gas arc) surrounding the central star 
cluster of NGC 604. 

2. OBSERVATIONS 

We conducted ^^CO(J=3-2) observations using the 
ASTE 10-m submillim eter telescope lo cated in the At- 
acama desert, Chile (jEzawa et al.ll2004[ ) from 2006 July 
to August. This was a part of our ADIoS project, which 
is an extragalactic CO(J=3-2) imaging survey to obtain 
a global view of the dense molecular medium in galaxies 
(iKohno et al.l[200l . We also conducted ^^CO{J^1~0) 
observations from 2005 December to 2006 March using 
the NRO 45-m telescope equipped with a 5 x 5 pixel focal- 
plane SIS array receiver (BEARS) capable of simultane- 
ously observing 25 positions in the sky (Sunad a et al. 
I2OOOI ) . Further information on the 45-m observations will 
be reported in detail in the forthcoming paper (Miura et 
al. 2007, in preparation). 

The front-end for the ASTE observations was a single- 
pixel cartridge-type 350 GHz SIS receiver (SC345; Kohno 
2005). An XF-type digital spectrometer was used to 
cover a velocity width of 445 km s~^ with a veloc- 
ity resolution of 5 km s^^ at 345 GHz. The observa- 
tions were conducted remotely from the ASTE oper- 
ation rooms in National Astronomical Observatory of 
Japan (NAOJ)-Mitaka and NRO, Japan, using the N- 
C0SM 0S3 network observati on system developed by 
NAOJ (|Kamazaki et al.l l2005l ) ^. The typical system 
temperature in double-sideband (DSB) was 200 K. The 
absolute pointing accuracy and main beam efficiency 
were verified by observing the CO(J=3-2) emission of 
o-Cet at interval of every 2 h. They were better than 2" 
r.m.s. and 0.6 during the observation runs, respectively. 
The stability of efficiency also was monitored using o- 
Cet, and was found to be stable within ±10%. 

The on-the-fly (OTF) mapping technique was em- 
ployed to obtain the ^^CO{J =3-2) and ^^CO{J =1-0) 
data. The "scanning noise" was removed by combin- 
ing the scan using the P LAIT algorithm described by 
lEmerson fc Graevel (I1988D . 

The full-width at half-power beams (FWHP) for the 
ASTE 10-m and NRO 45-m observations were 22" and 
16" at the rest frequencies of ^^CO(J=3-2) (345 GHz) 
and 12CO(J=1-0) (115 GHz), respectively We con- 
volved these maps to a common spatial resolution of 25"; 
this enabled us to measure CO ratios directly. 

3. RESULTS 

^ Observations using the ASTE were conducted remotely from 
Japan by using NTT's GEMnet2 and its partner R&E (Research 
and Education) networks, which are based on the AecessNova col- 
laboration of University of Chile, NTT laboratories, and National 
Astronomical Observatory of Japan. 



Fig. 1. — Three-color (B, V, and 6577Aband (Ha)) composite 
image of M 33 obtained using the 1.05-m Schmidt telescope in- 
stalled at KISO observatory. The image size is 59' X 55'. The 
white square represents the 5' X 5' area observed by the ASTE 
10-m and NRO 45-m telescopes; this area shows the location of 
the giant Hll region NGC 604. 

Figure 2 shows the velocity-integrated total intensity 
maps of 12CO(J=1-0) and ^^CO{J=3-2) emissions, in- 
tegrated over a velocity width of 200 km (from Vlsr 
= —350 km s^^ to —150 km s"^). Although this velocity 
range was larger than that of the emission, 80 km s~^, 
this has no effect on the total integrated intensity due to 
the the baseline subtraction. 

Wide-spread CO(J=1-0) emission can be seen within 
the mapped region. They are distributed around NGC 
604, on the arm to its north, and on the downstream 
side of the arm. The typical size of the clumps in the 
map is ~ 100 pc, and they aggregate to a larger complex 
that is similar to Giant Molecular Associations (GMAs; 

iRand fc Kulkarni I [1991 On applying the CLFIND 

method (IWilhams et al.l ll994D to the CO(J=1-0) data 
10 clumps were detected in this region and a few of them 
were new detections not reported in previous C0(J=1- 
0) observations. Detailed properties of these CO clumps 
will be reported in the forthcoming paper (Miura et al. 
2007). 

On the other hand, the ^^CO(J=3-2) map exhibits a 
compact morphology; the major CO(J=3-2) emission is 
confined to a small area close to the central cluster of 
NGC 604. Furthermore, the strong 1^00(7=3-2) peaks 
are located to the north of the i^CO(J=l-0) complex, 
i.e., the vicinity of the central star cluster of NGC 604. 
Interestingly, no clear CO(J=3-2) emission peak can 
be observed near the second-strongest C0( J=l-0) peak; 
this result reveals a significant variation in the physical 
properties of this region. 

The molecular clouds around NGC 604 are located to 
the south of the center of NGC 604. C0( J=3-2) emis- 
sion peak is closer to the central star cluster and sur- 
rounding the bright Hii regions. The southern side cor- 
responds to the upstream side in NGC 604, taking the 
galactic rotation of M 33 into consideration. 

We obtained the line ratio map of ^^CO{J—3~2) 
to i^CO(J=l-0) emissions, i?3_2/i-o- Figure 3 (left 
panel) shows the ratio of the integrated intensities of 
12C0( J=3-2) to i2C0( J=l-0). The ratio was computed 
for the position where the signal to noise ratios of both 
lines exceed 2. These ratios range from 0.3 to 1.3. The 
maximum it!3_2/i_o values in NGC 604 (1 - 1.3) are sim- 
ilar t o those in the center of the Milky Wav (lOka et al 



2007) and the central region of M 83 (|Muraoka et al.. 
2007t ). On the other hand, the lower values are simi- 



lar to those of the disk GMCs/GMAs in the Milky Way 
and the nearby sp iral galaxy M 31 (|Oka et al. 1 120071 : 
iTosaki et al.1 120071) . In fi gure 3, we found an arc-like 
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Fig. 2. — Velocity-integrated total intensity maps of the 
12CO(J=1-0) (left panel) and ^'^CO{J=3-2) (right panel) emis- 
sions in the 5' X 5' (1.3 X 1.3 kpc) region located in the center 
of NGC 604. The rms noise levels, a, in the CO(J=1-0) and 
CO(J=3— 2), were 0.5 and 1.3 K km s~^, respectively. The con- 
tour levels are 2, 4, 6, and 14o- for CO(J=1-0), and 1.5, 3.0, 
4.5, • • •, and 9a for CO( J=3-2). Note that the CO intensities are 
given in the main-beam temperature scale, i.e., J T^^,{v)dv. The 
dashed squares and crosses represent the region of the -R3_2/i_o 
map shown in figure 3 and the position of the central star cluster, 
respectively. These maps have a common spatial resolution of 25", 
indicated by a circle in the left-bottom corner. 
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Fig. 3.— Maps of the ratio of 1^00(7=3-2) to ^^CO{J=l-0) 
(left panel) and the Ha emission by HST (right panel; from the 
HST archive). The contour levels are 0.4, 0.6, 0.8, and 1.0. The 
crosshairs in the center indicates the position of the central star 
cluster of NGC 604, while the stars represent the peak positions 
of the H ll regions det ected by A 3.6-cm radio continuum emission 
llChurch wcll & Goss 1 993) . These maps clearly reveal the pres- 
ence of high ratio gas with an arc-like morphology surrounding the 
central star cluster of NGC 604. 

or shell-like distribution of a high-ratio gas surrounding 
the central star cluster of NGC 604. This "high-ratio 
gas arc" extends to the south-east to north-west direc- 
tion with a size of ~ 200 pc. 

4. DISCUSSION 
4.1. Nature of High-Ratio Gas Arc 

In this subsection, we discuss the physical conditions 
of the discovered high-ratio gas arc. 

A model calculation using large velocity gradient ap- 
proximation suggests that the high ratio gas {R3-2/1-0 ^ 
1) requires a kinetic temperature of Tkin > 60 K and a 
gas density of n(H2) ~ 10'^"'' cm"-^ (figure 9 in Muraoka 
et al. 2007). These values are in very good agreement 
with the results of Wilson et al. (1997) for the cloud 
NGC 604-2. These indicate that the warm and dense 
gases are distributed around the central star cluster with 
an arc-like morphology. 

In addi tion, there exists an arc-like d istribution of Ha 
emission (jGomez de Castro et al.l[200Clf ). wherein several 
compact Hii regions are em bedded, as revealed by ra - 
dio continuum observations (jChurchwell fc GossI Il999l ). 
These distributions exhibit a striking similarity to the 
western part of the high-ratio gas arc, suggesting that 
massive stars are formed within the warm and dense gas 
layers depicted as the high-ratio gas arc. 




Fig. 4.— Schematic view of NGC 604. 

Here it should be noted that the high i?3_2/i-o values 
in the arc are due to not only the high kinetic tempera- 
ture but also the high gas density. If the high -R3-2/1-0 
values are attributable solely to the UV heating from 
these newly formed stars, the high line ratio can only co- 
incide with the Ha shells, i.e., the high-ratio gas should 
be observed predominantly in the western part of Fig. 
3. However, our line ratio map indicates the presence of 
high-ratio gas, even in the eastern part of the molecular 
clouds, which is set apart from the bright Ha shell seen 
in Fig. 3 (right panel) . Furthermore, another calculation 
conducted using a photo dissociation region (PDR) sup- 
ports that the high fl3_2 /i-o values could be attributable 
to the high gas density (|Kaufman et al.iri999[ ). 

In conclusion, we suggest that the high-ratio gas arc 
reveals the presence of warm and dense molecular gas 
layers, and the massive stars now form within the dense 
gas arc. 

4.2. Triggered Dense Gas and Star Formation by 
"First-Generation" Star Formation 

In this subsection, we discuss the relationship among 
the high-ratio gas arc, the central star cluster, and the 
embedded young stars accompanying the Ha shell. 

The observational results obtained in this study can 
be explained by the following scenario. First, stars were 
formed in the northern part of GMA; this is referred to 
as "first-generation star formation" . These stars can now 
be observed as the central star cluster. Next, their stellar 
wind and supernova compressed the surrounding ISM. As 
a result, a dense gas layer was formed therein; such dense 
gases with high i?3_2/i_o ratios are distributed around 
the central star cluster. This is observed as the arc- like 
distribution of the high i?3_2/i_o ratio gas. New stars 
were then formed within the dense gases layer; this is re- 
ferred to as "second-generation star formation" triggered 
by the first-generation stars (i.e. the central star clus- 
ter). These secondary young stars formed in the dense 
gas layer can now be observed as the Ha shells and com- 
pact radio continuum sources as seen in Fig 3. Figure 4 
summarizes the schematic view of the proposed scenario 
for NGC 604. A similar situation is also reported in the 
Nil region, whic h is the second largest star forming 
region of LMC (Hatano et al.ll2006( l. 

Certain observations provide supporting evidence for 
the proposed scenario. First, the velocity field derived 
from the Ha emission of NGC 604 indicates that it 
consists of many filaments formed by multiple blowout 
events due to the energ y injected by massive stars 
(jTenorio-Tagle et al.ll2000l ). No Ha emission or CO emis- 
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sion peaks are seen in the central star cluster; however, 
the arc-like structure of the high i?3_2/i-o ratio gas has 
been associated with both Ha and CO emissions. These 
facts suggest that the Hii regions embedded in the high- 
ratio gas arc are in fact on- going star- forming regions, 
while the central star cluster that has already dispersed 
ISM is the "past" star-forming region. 

Second, the presence of shocks in the Ha shell was ob- 
serve d by optical spectroscopy (|G6mez de Castro et al.l 
I2OOOD : t hese are attributed to th e detected expanding 
motion (jTenorio-Tagle et al.ll200(i( ). This could be con- 
sistent with the view that the ISM surrounding the cen- 
tral star cluster was compressed by its energy blowout. 
The observation of a compac t Hi hole suggests the pres - 
ence of SNR in NGC 604 (|Deul k den Hartod fl99l : 
further, the observatio n of X-ray emitting hot plasma 
isanovic "eral][2006h may also support such expand- 
ing motion. 

Further evidence is provided by the optical 
narrow-band imaging of [Sii] and [Oiii] emissions 
(|Tenorio-Tagle et al.ll20M) . In the arc-like region with 
the high ratio, the [Olii]/Ha ratio is high while the 
[Sii] /Ha ratio is relatively low. Because both high 
[Oiii]/Ha and low [Sii]/Ha ratios correspond to a high 
excitation condition, it is suggested that the stars in 
the high-ratio gas arc region are significantly younger 
than those in other regions. This difference between 
the stellar ages of the central star cluster and the 
surrounding young stars is consistent with the proposed 
scenario. 

On the other hand, the age of the central cluster 
of N GC 604 has been estimated to be ^ 3 x 10^ 
years (jBruhweiler et al.ll2003t [Gonzalez Delgado &: Pered 
|2000[ ). Therefore, the timescale required to expand to the 
size of the detected Hii region, as indicated by the radio 
continuum observations, would be in the range of 10^ 
years to 8 x 10^ years (|Churchwell fc GossI [l999l ) . We 



found no significant difference between the central star 
cluster ages and the radio-derived ages. 

However, it should be noted that the stellar age was 
estimated over the region including both the high-ratio 
gas arc region and the central star cluster, i.e., the de- 
rived stellar age is a mixture of the proposed second- 
generation stars (in the Ha shell) and the first-generation 
ones (i.e., the central cluster). Consequently, inconsis- 
tency between these timescales and the proposed sce- 
nario is therefore not decisive and we need to conduct 
further observations to obtain spatially resolved stellar 
age measurements across the NGC 604 region. 

Note that the compressed warm dense gas is seen only 
in the southern and western parts of the star cluster. 
It is due to the direction from which material flows, 
namely, from the south-west to the north-east in spiral 
arm, and expanding material from first generation stars 
meets material flowing from the south-west into the arm. 
It suggests that a shock moves from the north-east to the 
south-west. 

Thus, we can conclude that the arc-like structure of 
the high-ratio gas is the site of second-generation star 
formation triggered by first-generation star formation. 
NGC 604 is an example of a large-scale sequential star 
formation. 
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